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A PREFACE TO THE NONADIABATIC THEORY OF 


ELECTRON-HYDROGEN' IONIZATION 

/ 

A. Temkin 

Laboratory for Theoretical Studies 
National Aeronautics and Space Administration 
Goddard Space Flight Ceiiter 
Greenbelt, Maryland 

Some preliminary considerations of the nonadiabatic theory of e-H 
ionization are given. The zeroth order problem is shown to give rise to 
an E^ threshold law. The stationary phase result for the asymptotic 
form of the zeroth order problem is shown to lead to a complete 
suppression of S-wave ionization events in which the energies of the 
scattered and ejected electron are equal. Extending the stationary 
phase results to our conjectured asymptotic form for the complete S-wave 
problem leads to an even greater suppression in the neighborhood of 
equal energy events. The combination of such an S-wave cross section 
together with reasonable higher partial wave cross sections yields a 
qualitative understanding of experimental ionization energy loss 
measurements in helium. 



I INTRODUCTION 


In this report, we continue to examine the very difficult problem of 
electron-hydrogen ionization from the point of view of the nonadiabatic 

theory'*'. In a previous note some preliminary results and observations were 

2 3 

stated (which were at variance with previous analyses of this problem ), 

which we shall attempt to quantify and amplify here. We shall continue to 

deal with total S-wave scattering. 

3 

Section II deals with the explicit demonstration of the E 2 ” threshold 

2 

law for the zeroth order problem . In section III we derive some of the 
salient physical consequences of the stationary phase integration for the 
closed form of the zeroth order ionization part of the wave function and 
extend the analysis to our conjectured asymptotic form for the complete 
S-wave problem. 

The final section discusses these results and how they may elucidate 
experimental results of the energy spectrum of electrons emerging from 
electron-helium ionization. Additional arguments are given for expecting 
a nonlinear threshold in the electron-atom ionization yield curve. 

For purposes of convenience, the main formulae of the nonadiabatic theory 
are given here. 
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'if is n r 2 times the S-wave function. Eq. (l.l) is the S-wave equation 
(energies in rydbergs, lengths in Bohr radii); Eq. (1.2) defines the basic 
expansion; Eqs. (1.3) are the coupled set of equations corresponding to (l.l). 
Eqs. (1.4) and (1,5) define the quantities occurring in (1.3). The symbol 
(Tmoo/no) is a Clebsch-Gordan coefficient. 
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II THE ZEROTH ORDER PROBLEM 


This problem is defined by neglecting the terms on the rhs of (1.3) for 
1=0. In the region ri ^ r 2 , this becomes: 


n' 2 ’ + dr 2 £ 


+ E 


) 


(o) 


(ri, r 2 ) = 0 


Expanding the solution in exact separable solutions, assuming E £ 0, yields 
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where 
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As E -» o from below this reduces to 
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E-«r 


= A sin Tx + 
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Consider the triplet case; the boundary condition for 


(o) 


is 


(ri, r 2 ) 


r x = r 2 = r 


(2.5) 
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Inserting Eq. (2.4) into (2.5) gives 


00 jL 

sin r Ris (r) = Y C e 11 R (r) 
L n ns 


+ f C(q)e _q - r F (r) dq 
J Q. 


The normalized radial functions are given by 


R ns 


) -4 e- r /“ F (l-»; 2J fO , 
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where F(a, b, c) is the confluent hypergeometric function. Expanding both 
sides of (2.6) is a power series in r and demanding that the equation be 
an identity in r, we find that the first two terms give rise to the equa- 
tions 


00 


l 


c 

n 

“ 3 - 



C(q) qr <iq 


o 


(2.8a) 


A = 




+ J C (q) q^ dq 

o 


(2.9) 


where (2.8) has been used to simplify (2-9). Let us now invoke the conser- 
vation of current (in the limit k -. — ) : 

n n 



( 2 . 10 ) 


Substituting (2-9) into (2.10) we arrive at 



- 8 - 



= 1 n l C n ’ 2 ( c ^)) ^ d( l (C(q))q^dq (2.1l) 

n=l 


If we assume that C n is of the form 


C = Ci/n^ , 
n ' 9 


( 2 . 12 ) 


then (2.1l) reduces to an identity when 


+ |i = 2^+1 ( ^ 


) 


providing 


(C(q)) q^ dq =xj£.(c 1 )J^ (C(q)) q^ dq 


(2.13) 
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Thus 


C 


n 



(2.14) 


satisfies the boundary condition to second order and is consistent with the 
conservation of current . 


The additional requirement (2.8a) provides no obstacle to C q 
being continuously extendable to ctfj) . If, for example, we assume 


COt) = - Ci e" Yq f 




then (2.8) is satisfied by a reasonable choice of y: 


V = [6/§(3)l* , 


where §(3) = 1.202057 is the Riemann zeta function. We have also 
examined the effect of satisfying the boundary condition to the 
next higher power in r. It turns out that all the previous equations 
can still be maintained and that one has one additional requirement 


to satisfy: 



Thus (2.12) still provides a consistent solution, but (2.15a) must 
be generalized to a two parameter extension in order that both (2.8a) 
and (2.8b) be satisfied. For instance the form 


cte) 


= -Cl e " Yl<3 - (q ? + b q^ ) 


can fulfill both equations. 

The fact that we have demonstrated that (2.14) provides a 
consistent solution is not a proof that it is correct. Eq. (2.1l) 
could be satisfied by not requiring the sum and integral parts be 
separately equal. We suspect, in fact, that the true solution is 
of this latter variety. The present demonstration is really in- 
tended to show that one can confidently expect (2.14) to be the first 
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term of an asymptotic series in inverse powers of n, and this is 

3 

all that is needed to complete the derivation of the E~ thr eshold 
law. 

The derivation of the singlet result is closely analogous to the 
above. The boundary condition is 




(o) 


Bn 


ri = r 2 s r 


0 , 


(2.16) 


! 

and the first two boundary condition equations are (2. 8k) plus 


CD 



J C(q)(q 2 ’ - q^) dq . 

o 


( 2 . 17 ) 



Substituting this into the conservation of current now yields 



(S - ?) dq 

o 


CO 



n=l 


Again this will be satisfied by (2.14) with the auxiliary condition that 
the q integral in ( 2 . 18 ) be zero. 

To show that (2.14) implies an threshold law we are indebted to 

L. 

K. Omidvar for a simple demonstration. Assume that one is considering 
excitation to a group of neighboring discrete (s) states for high n. The 
total cross section is given by 


n 2 

q(e) « y 


K l°J 


2 
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The energies of these states are given by 


e 


n 



( 2 . 20 ) 


from which it follows that there are 2n 3 d« (s) states in an energy range 

Converting (2.19) to an integral 


Q(E) « f k |C | 2 n 3 de , 

J n n ' n ' ( 2 . 11 ) 


continuing n into the continuum whereby 


k n = » 


( 2 . 22 ) 


and using (2. 14) for we obtain 


E 

Q(E) oc J vte~e n de n « iF , 

o 


(Z.Z2) 


Ill STATIONARY PHASE RESULTS AND THEIR PHYSICAL IMPLICATIONS 


Contrary to what was stated in Ref. 2 , the stationary phase result 
can be made to satisfy Kato's theorem, and therefore it may represent the 
correct asymptotic form of the zeroth order problem^. Explicitly for E > 0 
the separable solutions of (2-l) which do not vanish in the ionization region 
are ( + < 12 *^ = E ) : 


lim 

r l >r 2->°° 




e iqiri F 

12 


(r 2 ) dq2 


Stationary phase now gives for this integral 


tt (a) 

$ (°) _ fo(a) e i{\^3 p + — to (2 Je p)}^ 

ri>r 2 -.w ° p ^ 


lim 


wkev** 


W (*) - 
0 


I j /(U-'Vu 


oC 


r t >r^ 


\jc<&oc v , < r* . 

The requirement of Katols theorem (which is this case reduces to the con- 
tinuity of $ ^ and its first derivative at a = rr/4) is 
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to 


(a) 



9a 


(a) 



0 


( 3 - 3 ) 


Higher derivatives may have cusps. Likewise the assumption made for the 

2 

asymptotic form of the complete S-wave function leads under stationary 
phase to 


lim ^ = 
ri,r 2 -« 


i(q*.9ig) 



p + 


WcA«) ln(2v§ p) 


(3A) 


In the full S-wave case the validity of this result is very uncertain, 
hecause it requires a factorization of the exponential part through an in- 
finite sum of relative partial wave terms. However, if the expansion is not 
suitably convergent in relative partial waves or if the regions in which 
the asymptotic forms for the relative partial waves become valid differ 
sufficiently among themselves, it may be that such a factorization is not 
justified. Assuming it to be alright, then since Wo (a) has a cusp at 
a = for any value of 0 12 , whereas the true potential is completely 
continuous at a = J (®i2 i 0), it follows that the f(a ? © 12 ) and all its 



Nonzero functions 


-\ (d - 


derivatives (with respect to a) must be zero at a = 
of this kind can be constructed. For example. 


f(a.j0i 2 ) = g(©i 2 ) exp [- c/(l - -y)] 


(3-5) 


where 


r< rtan a 

^ " r> \cot a 


a 


< k 

. TT 

a > ij; 


(3.6a) 


The analyticity at a = | puts severe restrictions on the velocity 
distribution of the ionized particles. In the ionization cross section 
the quantity y is replaced by 


k 

< 


Y 


(3.6b) 


# 
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where k and k^ are the lesser and greater of k x and k 2 . It is convenient 

in exploring the physical consequences of (3*6) to work in terms a(e), 

the cross section for one of the electrons to appear with energy e . Such 

E 

a curve for a fixed available energy E must clearly he symmetric about — , 
thus 



(3-7) 


The quantity cr(e) is related to f by 


tU) = Iff 


ys 


E 

< 2 
E 

> 2 


(3-8) 


E 

In this way both a(e) and f are symmetric about 2 (and as long as f depends 

3 

on y and not E this gives rise to an iF threshold law) . Although the square 
root factor in (3-8) is essential to the threshold law, it has no important 
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effect on the shape of a as a function of r. In Fig. 1 we have suppressed 

TP 

the square root factor and have plotted {e £ ~ ) 



°/<Hr > 


( 3 - 9 ) 


for various choices of c. The normalization of a {e ) is chosen so that 

s 

% ( °) ■ 


IV DISCUSSION 
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The primary observation to be made from Fig. 1 is that the region of 
momentum space in which the two outgoing particles have equal speeds is 
highly suppressed. This is an obvious consequence of the zeroth order prob- 
lem in which any correlations which could be otherwise taken up in O12 
are eliminated in the definition of the model (thfis all correlations are 
necessarily in relative energies only). It is an interesting consequence 
of the full S-wave problem with the asymptotic form of Eq. (3-4) that this 
suppression is enhanced. 1 Such an unexpected consequence might in fact be 
used to argue against Eq. (3-4) (.and we do not think that (3*4) can rigour- 
ously be correct) and possibly against the asymptotic form (15) of Ref. 2. 
However, in the absence of a definitive analysis, it is worthwhile to be 
guided as much as much as possible by experiment. 

Before doing this, let us note that distributions of Fig. 1 are 

directly opposite to what one would expect from phase space. For phase 

1 £ 

space is proportional to [e(E-e)] 2 which has its maximum at e = ~ and is 
zero at the end points. 

In Fig. 2 we have plotted the energy loss cross section, a(e), for S-wave 
together with assumed P and D waves cross sections. The experiment that is 
envisioned here is one in which one of the emerging particles in the ioniza- 
tion process is observed with energy e in the forward direction relative to 
an incoming beam of energy k 2 = E + 1^; 
electrons in the atom. The energy of the second particle is therefore E - e, 
and the experiment automatically averages over the directions of the second 
particle relative to the first (©is). The S-wave is independent of the 


1^ being the binding energy of ejected 
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angle of observation relative to the direction of the incoming beam, there - 

E 

fore it is necessarily symmetric in e about — for any angle of observa- 
tion. This is not so for higher partial waves in which cases only the 
integrated cross sections over all angles of observation is symmetric in 
e. Appendix I contains a brief analysis of the P wave. For zero angle of 
observation the higher partial cross sections are drawn as increasing func- 
tion of e corresponding to the simple physical expectation that the more 
energetic the collision products, the more in the forward direction they 
should emerge. In the Appendix it is shown why the p wave, however, must 
decrease for e near its maximum value of E. The curve representing the 

sum of these partial cross sections is in qualitative accord with preliminary 
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experimental results of He i deman , for the case of electron-helium ionization; 
The above experiment is not at low enough energies, as compared with 

g 

that of McGowan et al. , to establish the form of the threshold energy 
dependence of the ionization cross section. Rather it serves to illustrate 
the primary correlation of the two outgoing particles as an anticorrelation in 
the available energy. However, this correlation itself can be quite signifi- 
cant in determining the form of the threshold law. To illustrate this 
point, let us note that when we say that in elastic scattering the scattered 
electron is completely shielded from the nucleus we do not mean 


lim 

rx-«» 


Y = sin (krx + 6) Rx s (r 2 ) 


( 4 . 1 ) 



r 


~ 2 f - 


for all values of r g . In particular it is probably incorrect as r 2 but 

that makes no difference, since the wave function itself vanishes in the joint 
limit. In the case of ionization, the region ri a* r 2 -» 00 is not energetically 
disallowed, however, it may be dynamically unfavored to the extent of pro- 
viding complete shielding for the outer particle. That is what, in effect, 
we have been arguing happens. 

In addition to the dynamical effects that we have already discussed, 

there are additional effects which we believe should strengthen the cogency 

of the concept of shielding. It has been customary in dealing with this 

problem to divide the asymptotic region into a reaction zone and an emerging 
o 

zone or zones . This implies that the phenomenon or ionization takes place 
fairly close to the nucleus and that as the particles emerge their "orbits" 
change relatively slowly as a result of the asymptotic interactions. However, 
at any distance from the nucleus there are bound orbits, and it always is 
possible that the inner particle, in particular, as it emerges can get 
caught in a bo und state. If such transitions can occur into high bound 
states, they can also occur to states of lower continuum energy also. Thus 
the combination of possibilities can only serve to separate further the 
energies of the emergent particles, and to lower the absolute ionization 
cross section as a function of energy relative to elastic and inelastic 
events. Near threshold this is an additional argument for a nonlinear 
(concave up) dependence. 



APPENDIX 


We give here a brief outline of , the analysis of P-wave ionization. 

Let 9-n, be the spherical angles of an observed electron in an ionization 

event relative to a fixed coordinate system whose z-axis in the direction 

10 

of the incoming electron. The P-wave function can be written 


Y p = cos 0i2 [f p + cos 0 12 fpl 


+ sin 9 fi cos y b sin 0 12 f p } 


where for purposes of computing the ionization amplitude 
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Y-n is angle between ki~k 2 plane and the z-axis. 

-b 

Actually the fp should be multiplied by a function of ri and r 2 which 
gives outoing current for two particles in the directions ^2 respectively. 
Although there is some discussion as to what this function should be, wc 
shall assume that for higher partial waves, it is not such as to exclude 
events in which ki = k 2 . 

The probability for a given event characteristed by the two vectors ki 
and is the absolute square of Yp- For one electron to be observed in the 
direction 0^ and the second to be anywhere, we have 


<Tp(e & ) = J I Yp(kx k 2 © B sin 9 12 d 0 12 dy B 


Special cases: 


0 B = 0 


TT 

Op(°) = 2 I (t f pl + l f p! cos2 0 12 J sin 0 12 d0 12 
o° 
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j % [fpfp*] Sin 0 12 d0 12 
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B 


TT 


2 


^ 2 _ pTT ~ 2 

°p (| y = i j sin3 012 d012 

o 

Note also that <jp(0) = op(Tr). Thus if it is unlikely to find a high 

e (p - wave} particle in the backward direction, then it must also be unlikely 

to find one in the forward direction. This is the justification for having 

the P-wave descend rapidly for the higher e as pictured in Fig. 2. 

Although the various cr-n(©-n ) cross sections for a given ©„ are not 

r n Jd 

E 

symmetric about e = g , we show below that the integrated cross section 

over all Q_ is symmetric, as it must be, since each ionization event which 
Jd 

gives rise to an electron of energy e must also give rise to one of 
energy E-e somewhere on the scattering sphere. 

Integrating over the whole sphere, we find 


- - 



sin Q_ 

15 


d8 B 
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sin 0 12 d© 12 


fClTpl + I f p I 3 ] sin ©12 


d© 


12 


j* 


(f f *) 
P ' 


sin 2© 12 &©i 2 


The first integral is, recalling (A 2 ), obviously symmetric with respect to 
ki 2 k 2 . If one separates f p and f p into real and imaginary parts and uses 
the symmetry under ki k 2 of each part separately, one can see that the 
second integral is also symmetric. 

Finally for higher partial waves a similar analysis can be carried out 
by using symmetric Euler angles 1 '*' and converting to the Hylleraas -Breit 
Euler angles (labelled by the subscript B above) via Eqs. (3.l)and (3.2) 
of Ref. H, 
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FIGURE CAPTIONS 


S-wave ionization cross sections vs. the residual energy e for 
a fixed bombarding energy calculated on the basis of Eq. ( 3 * 9 )> 
for various values of c. The experiment of Heideman (Ref. 7 ) 
indicates a large value of c. 

S-wave cross section corresponding to c = 2 of Fig. 1. The 
higher partial waves are supposed to correspond to ionization 
events in which one of the emerging particles is observed in the 


forward direction. 
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